Field induced phase transitions in the helimagnet Ba 2 CuGe 2 07 



(N 

o 

> 
O 

in 

(N 



u 

-f— > 

-4-* 
03 



i 

c 

o 
o 



> 

(N 
00 

in 



(N 



J. Chovan 1 '*, M. Marder 2 ^ and N. Papanicolaou 3, * 
1 Department of Physics, Matej Bel University, Banska Bystrica, Slovaki^\ 
2 Center for Nonlinear Dynamics and Department of Physics, 
The University of Texas at Austin, Austin, Texas 78712 Q and 
3 Department of Physics and Institute of Plasma Physics, University of Crete, Heraklion, Greec^\ 

(Dated: November 27, 2012) 

Abstract. We present a theoretical study of magnetic properties of the two-dimensional spiral 
antiferromagnet Ba2CuGe2 07 in the presence of an external magnetic field. Some time ago we pre- 
dicted that when a field is applied parallel to the c axis, a previously anticipated Dzyaloshinskii-type 
incommensurate-to-commensurate phase transition is actually mediated by an intermediate phase. 
This phase was recently observed in experiments where magnetic fields were nearly parallel to the 
c axis and referred to as a double-k structure. Here we calculate the low-energy magnon spectrum 
throughout the phase transition, verifying stability of the intermediate phase, and providing the 
opportunity for comparison with inelastic neutron scattering. In addition, we calculate the T — 
phase diagram for arbitrary canted magnetic fields, in general agreement with experiment. In partic- 
ular we show that a weakly broken U(l) symmetry of our model accounts for sudden 7r/2 rotations 
of the magnetic structures observed in experiment. Finally, our analysis suggests a nonzero weak- 
ferromagnetic component in the underlying Dzyaloshinskii-Moriya anisotropy, which is not crucial 
for the appearance of the double-fc structure, but is important to restore quantitative agreement 
with experiment. 

PACS numbers: 75.30.Ds, 75.30.Gw, 75.30.Kz 



I. INTRODUCTION 

The presence of Dzyaloshinskii-Moriya (DM) 
anisotropy 1,2 in low-symmetry magnetic crystals 
typically leads to weak ferromagnetism, as a result of 
slight spin canting in an otherwise antiferromagnetic 
(AF) ground state. Another possibility is the occurrence 
of helimagnetism whereby spins are arrayed in a helical 
or spiral structure whose period (pitch) extends over 
several decades of unit cells. Ba2CuGe2C>7 is an example 
of a helimagnet which has proved to be especially 
well suited for experimental investigation thanks to a 
fortunate combination of physical properties. It is an 
insulator whose magnetic properties can be understood 
in terms of localized s = i spins carried by the Cu 2+ 
ions. The scale of energy set by an exchange constant 
J 1 meV is very convenient for neutron scattering 
experiments. Because of the low tetragonal symmetry 
(space group P42im) the corresponding Heisenberg 
Hamiltonian involves an interesting combination of 
antisymmetric (DM) as well as symmetric exchange 
anisotropics which lead to a rich phase diagram. In 
particular, the strength of anisotropy is such that 
magnetic phase transitions take place at critical fields 
that are well within experimental reach. 

Indeed, a series of experiments in the late nineties 3 - - — 
revealed the existence of a Dzyaloshinskii-type 8 
incommensurate-to-commensurate (IC) phase transition 
when the strength of an external field applied along the 
c axis exceeds a critical value, H c ^2 T. For H < H c 
the ground state is an incommensurate spiral whose 
period L = L(H) grows to infinity in the limit H — > H c . 
For H > H c the ground state was thought to become a 
commensurate antiferromagnetic spin-flop state. We also 



note that the Dzyaloshinskii-type transition is similar to 
the cholesteric-nematic phase transition induced by an 
external magnetic field in chiral liquid crystals^—. 

However, a detailed theoretical investigation^ 2 - - — in- 
spired by the earlier experimental work 3 - - — predicted that 
the IC phase transition does not occur immediately, and 
that, instead, between the incommensurate and commen- 
surate phases occurs a separate intermediate phase. In 
short, there exist two critical fields, H c \ and H C 2 such 
that H c i < H c < H C 2 where H c ~ 2 T is the critical 
field for the presumed Dzyaloshinskii-type phase transi- 
tion. For H < H c i ~ 1.7 T the ground state is a flat 
spiral (cycloid) that propagates along the x axis while 
the staggered magnetization rotates in the xz plane. For 
H > H c i the cycloid transforms into a nonflat spiral 
where all three components of the staggered magnetiza- 
tion are different from zero (Fig. [I]). Such a state may 
concisely be described as an antiferromagnetic conical 
spiral that propagates along the x axis while it nutates 
around the y axis. When the field approaches a second 
critical value, H c2 ~ 2.9 T the antiferromagnetic cone 
degenerates into a uniform state whose staggered mag- 
netization points along the y axis modulo an azimuthal 
rotation. This state is a commensurate antiferromagnetic 
spin- flop state which is the ground state for all H > H C 2- 
Therefore the Dzyaloshinskii field H c is not a true crit- 
ical field, and the corresponding IC phase transition is 
actually mediated by an additional phase in the region 
H c i < H < H C 2- 

This prediction 12 remained unexplored for almost a 
decade. However a new series of experiments has now 
confirmed the occurrence of an intermediate phase in 
the form of an antiferromagnetic conical spiral which has 
been called a double-k structure by the experimental dis- 
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Figure 1. (Color online) Diagrams illustrating the forms of 
various spiral states for differing values of magnetic field h 
in units of f .68 T along the c axis. Blue lines on the sphere 
surface trace out directions for the staggered magnetization, 
placing the base of each n at the center of the sphere and then 
moving from one unit cell to the next along x. the direction 
of spiral propagation. For < h < h c \ the spins describe a 
cycle in the xz plane, while for 1 < h < h C 2 = v3 they also 
have a nonzero oscillating y component. 



coverers. This state occurs as predicted when an external 
magnetic field is applied almost perfectly parallel to the 
c axisi^, while further experiments have also explored 
the T = phase diagram in the presence of an arbitrary 
canted magnetic fielcU^. Our current task is to confirm 
that the recently observed double- k structure is indeed 
the intermediate phase predicted in Ref. 12 and further 
to calculate the phase diagram in arbitrary canted mag- 
netic fields so as to complete the connection with the 
latest experiments. 

In the method of calculation we closely follow the work 
of Ref. 12. We first note that Ba2CuGe207 is a layered 
compound where the Cu atoms form a perfect square 
lattice within each layer. The interaction between out- 
of-plane neighbors is ferromagnetic and weak. Therefore 
the inter layer coupling is ignored in the following dis- 
cussion which concentrates on the two-dimensional dy- 
namics within each layer. The low-energy dynamics is 
then calculated in terms of a suitable variation of the 
nonlinear a model governed by Lagrangian density C : 



C = C - V- (1) 

Co = 2^° n ' + h ■ n x don; 

1 2 1 2 

V = - [din - e 2 x n) + - (d 2 n -e l xn) 

+ i(n-h) 2 +d 2 (hxe 3 )-n. 

Here ei, e 2 , and e 3 are unit vectors along the x, y, and 
z axes and n = niei + n 2 e 2 + ^3^3 is the staggered mag- 



netization which is a unit vector field that depends upon 
the in-plane spatial coordinates x and y as well as the 
time variable t: n = n(x,y,t). Accordingly, derivatives 
are described by d\ = d/dx, d 2 — d/dy, and do = d/dt. 
The applied magnetic field h = h\e\ + h 2 e 2 + ^3^3 may 
point in any arbitrary direction. It should be noted that 
anisotropy has been restricted to the special KSEA limit 
by setting n = in the Lagrangian of Ref. 12. Oth- 
erwise Eq. ([T]) gives the most general Lagrangian com- 
patible with symmetry, expressed in fully rationalized 
units. Specifically, frequency is measured in units of 
tko = 0.24 meV, distance in units of 33.75 A and mag- 
netic field in units of 1.68 T. The remaining free param- 
eter d z — \/2D z /D± is a dimensionless ratio of the out- 
of-plane component of the DM anisotropy D z and its 
in-plane component Z?^.The latter has been completely 
suppressed in Eq. ([T]) except through the definitions of 
rationalized units quoted above. 

The ground state properties and the associated low- 
energy dynamics will be calculated from Eq. (|TJ) for a 
magnetic field of varying strength and direction. Hence, 
in Sec. |TT] the field is restricted to point along the c (or 
z) axis and its strength is varied through the IC transi- 
tion. We recover the results of Ref. 12 and complete an 
explicit calculation of the low-energy magnon spectrum 
throughout the intermediate phase. In Sec. IHII we study 
the case of a field applied in a direction perpendicular to 
the c axis. We thus recover the experimentally observed 
bisection rule and further illuminate the role of the out- 
of-plane DM anisotropy d z . The case of a magnetic field 
applied in an arbitrary direction (canted magnetic field) 
is analyzed in Sec. II VI where we present a theoretical pre- 
diction for the T = phase diagram, in fair agreement 
with recent experiments. Finally our main conclusions 
are summarized in Sec. [V] 



II. FIELD PARALLEL TO c 

We begin by specializing to the case where the mag- 
netic field is applied strictly along the c axis: h = he^. 
Then the potential V of Eq. (JT|) reduces to 



(dmy + (d 2 ny + (1 + h 2 ) n\ + 1 
- [(dxni - d 2 n 2 ) n 3 - (nidi - n 2 d 2 ) n 3 ] (2) 
and is symmetric under the U(l) transformation 

x + iy — >• (x + iy) e 1 ^ , rii + in 2 (m + in 2 ) e~ , (3) 

which is somewhat unusual in that an azimuthal rotation 
of spatial coordinates x and y by an angle ip is followed by 
a corresponding rotation of the staggered magnetization 
by an angle —ip. 

The ground state is obtained by finding energy- 
minimizing solutions n of the static energy functional 



/ 



W = / dxdyV. 



(4) 
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Figure 2. Solutions of Eqs. (0 minimizing the average energy density of Eq.fTS]) for a number of illustrative magnetic fields 
h — h z pointing along the c axis. Note that the period L varies with the field. 



In order to enforce the constraint that the staggered mag- 
netization be of unit length we adopt a parameterization 



n = sin <f> sin O ei + cos O e2 + cos <f> sin e3, (5) 

which differs from more standard parameterizations by a 
circular permutation, but turns out to yield slightly more 
compact expressions later on. 

Our first task to to find solutions = Q(x,y), $ = 
<fr(x,y) that minimize W. Extensive two-dimensional 
simulations yielded negative results for a potential 
ground state in the form of, say, a vortex lattice. On the 
other hand, nontrivial results are obtained by restricting 
both and $ to be functions of x alone: 



S{x,y) =6{x); 9(x,y) = </>(x). 



(0) 



In view of the U(l) symmetry in Eq. Q any solution 
we find of this type automatically produces a family of 
additional solutions of the same energy rotated by angle 
ip. Varying W then yields 



oh 



(2 di<j> — 2) cos 9 d\0 + j 2 cos cj> sin</> sin 9 



sm( 



die = ((<9i</>) 2 - 2d l <p + 7 2 cos 2 ^ cos (9 sin0. (7) 

with 7 2 = 1 + h 2 . Here subscript 1 indicates a derivative 
with respect to x. All derivatives with respect to y vanish 
because we are working in a space of one-dimensional 
solutions. 

To illustrate the solutions, we first consider the special 
case of a flat spiral (cycloid) with 9 = tt/2. Then the 



second of Eqs. ([7|) is automatically satisfied and the first 
becomes 

<9 2 + 7 2 cos0sin</> = 0, (8) 

while the staggered magnetization becomes 

n = (sin</>, 0, cos^>), (9) 

a cycloid that propagates along the x axis while rotating 
in the xz plane (upper left panel of Fig. [l]). The solution 
for 4> obeys 



di4> = \/ ft 2 + 7 2 cos 2 0, x = 



o \/ S 2 + 7 2 cos 2 ip 



(10) 

The result can be expressed in terms of elliptic functions 
but there is no particular advantage to doing so. S 2 is 
a positive constant that will be determined below. The 
cycloid has a period (pitch) of 



L = 



2tv 



o \/S 2 + 7 2 



(11) 



and the free parameter S is determined by the require- 
ment that the average energy density w — W/ L achieve 
a minimum: 

— / d<t>^8 2 + 7 2 cos 2 4>=l^w = -(l- S 2 ). (12) 
2tt J q 2 

As 7 (or h) increases S becomes zero at a critical field: 
7 = 7c = tt/2 =>h = h c - 



lwl.21. (13) 
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Figure 3. Magnon spectrum along the x direction in the re- 
duced zone scheme for six illustrative magnetic fields. For 
h > h c i = 1.01 the results show the magnon spectrum of the 
intermediate phase. The wave number Q\ is measured in rela- 
tive lattice units defined in Ref. 12. Note that the lowest-lying 
band has a linear dispersion relation. 

In physical units, H c = 2.04 T. This is the Dzyaloshin- 
skii critical field and the corresponding Dzyaloshinskii 
scenario may be described as follows: for h < h c the so- 
lution is a flat spiral that propagates along the x axis and 
rotates in the xz plane. As h approaches h c the spiral is 
highly distorted and becomes a kink-like structure with 
diverging period. For h > h c the ground state becomes 
the uniform spin-flop state 

n= (1,0,0) modulo 17(1). 

We realized that this scenario was incomplete when we 
computed the magnon spectrum^ of the flat spiral and 
found negative eigenvalues starting at 

h cl = 1.01, Hci = 1.7 T. (14) 

Above this value the flat spiral is unstable. We thus 
return to energy minimization and revoke the assump- 
tion 9 = 7r/2, although continuing to assume a one- 
dimensional structure of the form cj> = 4>{x) and 9 = 9{x). 

Efforts to find explicit analytical solutions of Eqs. ([7]) 
have not been fruitful so we resort to numerics. We min- 
imize the energy density 

w = ~ [ dxV{0,4>) (15) 
L Jo 
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Q 1 [r.l.u.] Q 1 [r.l.u.] Q 1 [r.l.u.] 

Figure 4. Magnon spectrum along the x direction in the ex- 
tended zone scheme for six illustrative magnetic fields. For 
h > h c i = 1.01 the results show the magnon spectrum of the 
intermediate phase. Bands have been assembled in a fashion 
that corresponds with conventions in publication of experi- 
ments. 



over a periodic chain of length L and vary L to achieve 
a minimum for any given value of h. 

For h < h c i — 1-01 we recover the previous results for 
the fiat spiral. But for h > h c \ a nonflat spiral arises 
with nontrivial <fr(x) as well as 9(x). We emphasize again 
that our numerical investigations of more general two- 
dimensional structures have never yielded lower energies. 
Examples appear in Fig. [5] for a variety of field values. 
Entering the intermediate phase for h > h c \ , ri2 acquires 
nonzero values and one can describe the state as an anti- 
ferromagnetic conical spiral that propagates along x but 
nutates around y. This is precisely the structure de- 
duced from recent scattering experiment a 16 i 17 and called 
a double-k structure because of two-fold peak character- 
istically observed during experimental scans through k 
space. As h increases, the component 712 becomes larger 
and larger until at h C 2 = v3 (or H C 2 — 2.9 T) the so- 
lution becomes a spin- flop state with n = (0, 1, 0). This 
upper critical point was determined in Ref. 12 from a sta- 
bility analysis of the spin-flop state. The existence of the 
intermediate state does not depend upon the presence of 
a nonzero transverse magnetic field. 

We now proceed to obtain the magnon spectrum of 
the intermediate state. We consider perturbations in the 
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Figure 5. Magnon spectrum along the y direction for six 
illustrative magnetic fields. For h > h c \ = 1.01 the results 
show the magnon spectrum of the intermediate phase. The 
wave number Q2 is measured in relative lattice units defined 
in Ref. 12. 

form 

®(x,V,t) = 9{x) - g(x,y,t) 

$(x,y,t) =4>(x)+f{x,y,t)/sin9(x) (16) 

where 9 and <f> are solutions for the intermediate state 
found previously. Then working to linear order in /, g 

(df + d\ - d 2 )f = U u f + U 12 g + Ad l9 + Bd 2 g + Cd Q g 
(Of + dl - dl)g = U 22 g + U 21 f - Adxf - Bd 2 f - Cd a f 

(17) 

where all functions except / and g are functions only of 
x and are given by 

U n = - (di9f + cos 2 9 ({d^f - 2 di<fy 
+7 2 (cos 2 <f> cos 2 9-2 cos 2 <j> + l) 

{2d 1 <t>-2)d 1 9 

Ul2 — :— 7 

smtf 

(2 d\4> — 2) cos 2 9 d\9 + 2~f 2 cos</> sin</f cos9 s\n.9 
U 2 i = — 

U 22 = - ((c)^) 2 -2d 1( f> + 7 2 cos 2 (2 sin 2 9-1) 

A = (2d 1 (t>-2) cos 9 



B = — 2 sin (j) sin 9 

C = -2h cos <j) sin 9. (18) 

We have verified that for the flat spiral (9 — n/2, 
d\4> — S 2 + 7 2 cos 2 4>) these expressions reduce to those 
previously obtained for the magnon spectrum of the flat 
spiral in Eq. (5.2) of Ref. 12, but with </» i — > 9. 

We have solved the linear system (|17p by a Bloch anal- 
ysis of the type given in Appendix A of Ref. 12 now 
extended to calculate the low-energy magnon spectrum 
throughout the intermediate phase h c % < h < h c2 . We 
present the results of the magnon calculations in Figs. [3] 
through [5] We make the following comments: 

• All eigenvalues are positive. Therefore the inter- 
mediate state is locally stable. This computation 
does not prove it is the ground state, but in com- 
bination with extensive numerical explorations of 
two-dimensional states that found no solutions of 
lower energy, it is a strong indication. 

• We provide plots both in the reduced zone scheme 
and the extended zone scheme. The reduced zone 
scheme is more compact, particularly for h c \ and 
below. However as the field increases towards h c2 
the reduced zone scheme acquires a large number 
of bands that are resolved more clearly in the ex- 
tended zone scheme. Experimentalists are likely to 
find the display in the extended zone scheme more 
useful. 

• Along Q\ the low-energy spectrum is linear at the 
zone center. Moving towards h c2 it acquires two 
bands, an 'acoustic' band with linear dispersion and 
an upper optical band (higher bands exist that have 
not been resolved by the computation). The linear 
portion of the acoustic band is the Goldstone mode 
of these magnetic spin states. In the limit that 
h — s- h c2 the bands depicted here collapse onto the 
horizontal axis; the next excitation is at an energy 
over 0.4 that lies above the top of the figure. 

• Along Q 2 the low-energy spectrum is quadratic. As 
h increases towards h c2 the quadratic regions be- 
come small and the spectrum becomes nearly lin- 
ear. Upon reaching h c2 , the dispersion becomes 
completely linear. At this point it produces the 
Goldstone mode of the spin-flop phase. 



III. FIELD PERPENDICULAR TO c 

We next consider a field applied in a direction strictly 
perpendicular to the c axis, a case that had attracted 
experimental interest already in Ref. 6. For the moment, 
we assume that the field is applied along the y axis, h = 



6 



(0, /ij.,0), hence the potential V of Eq. ([T]) reduces to 



V 



(chn) z + {d 2 ny + n 2 3 + hinj + 1 



+ h±_d z n\ 



[{dini - d 2 n 2 ) n 3 - (nidi - n 2 d 2 ) n 3 ] 



(19) 



where the applied field enters in two distinct ways; 
namely through the appearance of an effective easy-plane 
anisotropy \h\ri 2 and a Zeeman-like anisotropy h\_d z ii\. 
The latter also contains the strength d z of the out-of- 
plane oscillating component (±D Z ) of the DM vectors 
which was neglected in the analysis of Ref. 6. 

To find minima of the energy functional we first note 
that the positive term ^h\n\ again favors a flat-spiral 
configuration with n 2 — which propagates along the x 
axis. Using the angular parametrization we write 



which is inserted in Eq. (|19l) to yield 

{di<t>- if +cos 2 
where the only free parameter 



h sin ( 



(21) 



h i d. 



(22) 



is a combination of the applied field h± and the effective 
out-of-plane DM anisotropy d z . 

Otherwise, the calculation is similar to that of the flat 
spiral in Sec. II. Stationary points of the energy func- 
tional W = j Vdx now satisfy the ordinary differential 
equation 

d 2 <fi + cos (p sin — h cos (j) — (23) 

whose first integral is given by 

(dicj)) 2 - cos 2 - 2fein0 = C = 2h + S 2 . (24) 

Our choice of the integration constant C indicates that 
minimum energy is achieved with a positive new constant 
denoted by 8 2 . The actual configuration $ = (j> (x) is then 
given by the implicit equation 



dip 

^S 2 + cosV + 2/i(l + sinip) ' 
and the corresponding spiral period L is given by 

L -- 



(25) 



(26) 



yj5 2 + cos 2 4> + 2h{\ + sindi) 

Finally, the free parameter 5 2 is calculated by minimiz- 
ing the average energy density w — ^ Jq V(x)dx which 
yields 



— / d(j>JS 2 + cos 2 + 2h(l + sm<f>) = 1 
2n In 



(27) 




Figure 6. Illustration of the bisection rule in a transverse 
field applied along the y" axis. The spiral propagation vector 
points along the x' axis, while the staggered magnetization 
rotates in the x" z plane. 



an algebraic equation that may be used to determine S 2 
for each value of h. The corresponding minimum energy 
is then given by 



1 



w 



(1-5 2 - 2h). 



(28) 



In the absence of the out-of-plane DM anisotropy 
(d z = 0) the configuration just calculated reduces to the 
zero-field flat spiral of Sec. II for any value of the applied 
transverse field because h = h±d z — for all h±. In par- 
ticular, no phase transition of the Dzyaloshinskii type 
would be expected to occur for a field applied in a direc- 
tion strictly perpendicular to the c axis, as presumed in 
the analysis of early experiments^. 

However, the situation changes significantly for d z =^ 
0. Then the effective field h = h±d z is different from 
zero except when h± — 0. With increasing hx, and thus 
increasing h, the parameter S 2 decreases and eventually 
vanishes when h reaches a critical value h — h c computed 
from Eq. (|2"6l applied for S 2 — 0. A simple numerical 
calculation yields h c = h c 1 d z = 0.3161, or 



0.3161 



(29) 



In the limit h± —> h\, 8 vanishes and the average en- 
ergy density of Eq. (|2"8"|) reduces to io = |(1 — 2h c ) 
which coincides with the energy of the uniform spin- 
flop state n = (—1,0,0). Thus we again encounter a 
Dzyaloshinskii-type phase transition at a critical field 
that now depends on d z . 

As mentioned already, no such transition was detected 
in the early experiments 6 which were conducted with 
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transverse magnetic fields of limited strength H± < 2 
T or h± < 2/f .68 « f .2. However, recent experiments 17 
reveal a critical field = 9 T or h\ = 9/1 .68 = 5.36 
and, using Eq. (|2"9"|) . 



0.06. 



(30) 



As far as we know, this is the first estimate of the strength 
of the out-of-plane DM anisotropy and will be used in 
all numerical calculations presented in the continuation 
of this paper. Incidentally, using the definition of the 
rationalized anisotropy d z = \/2D z /eJ from Ref. 12, we 
find D z /J = 0.0076, to be compared with e = D±/J = 
0.18. 

The preceding calculation was completed in Ref. 13 
with a detailed calculation of the corresponding magnon 
spectrum which could prove useful for the analysis of fu- 
ture inelastic neutron scattering experiments in the pres- 
ence of a strong transverse magnetic field H± . The same 
calculation reveals no sign of further critical instabilities 
as long as d z < 0.5. In particular, an intermediate phase 
of the type encountered in Sec. II is not present in the 
case of strictly transverse magnetic fields and d z < 0.5. 
In other words, the predicted phase transition is of pure 
Dzyaloshinskii typeS. 

This section is completed with a brief discussion of the 
case of a transverse magnetic field 



= h± (sin ip, cos ip, 0) 



(31) 



which points in an arbitrary direction within the basal 
plane obtained by a clockwise rotation of the y axis 
with angle ip (see Fig. 6). In fact, the ground-state 
configuration for this more general case (ip ^ 0) can 
be surmised from the special ip = solution calcu- 
lated earlier in this section by simple algebraic trans- 
formations, thanks to the underlying U(l) symmetry of 
Eq. © broken by the applied transverse field. Indeed, 
let rii — ni(x),n 2 = 0, 713 = 77,3(2;) be the ip = solution. 
Then the solution for ip ^ is given by 

n[ = cosipm(0 , n' 2 = -sinV>ni(£) , n' 3 = n 3 (£) , 

where £ = xcostf) + ysiaip. Thus the new spiral propa- 
gates along the x' axis obtained by a counter-clockwise 
rotation of the x axis with angle ip (see Fig. 6) while the 
staggered magnetization rotates in the plane x" z which 
is perpendicular to the field direction (axis y"). In other 
words, a flat spiral (cycloid) that initially propagates 
along the x axis and rotates in the xz plane (tp = 0) 
is reoriented to propagate along the x' axis (ip 7^ 0) so 
that the normal to the spin plane (axis y") points along 
the applied magnetic field. The angle formed by the di- 
rection of spiral propagation (axis x') and the normal to 
the spin plane (axis y") is bisected by the conventional 
crystal axis b — (0, 1, 0) denoted by a dotted line in Fig. 6 
for any ip. When the field is applied along b, tp = ^ and 
the normal to the spin-rotation plane is parallel to the 
propagation vector (screw-type spiral) . 



The "bisection rule" just described theoretically was 
experimentally discovered already in Ref. 6. Actually, 
agreement with the ideal bisection rule requires that 
H±_ > 0.5 T in order to overcome a certain energy 
barrier due to discreteness effects which lead to an ad- 
ditional tetragonal anisotropy that breaks the underly- 
ing U(l) symmetry even in the absence of a transverse 
field&i£. The same anisotropy explains the experimen- 
tal fact that the spiral propagates along the x = (1,1,0) 
or x = (1,1,0) directions, in the absence of a transverse 
field, while a sufficiently strong field H±_ > 0.5 T is re- 
quired to reorient the spiral according to the bisection 
rule. 

We have thus completed the discussion of the phase 
diagram in the presence of a field strictly parallel to 
the c axis (Sec. II) or a field strictly perpendicular to 
c (Sec. III). The general case of a canted magnetic field 
is discussed in the following Sec. IV. 



IV. CANTED MAGNETIC FIELDS 

We now turn our attention to the most general case of 
the applied field h, whose transverse component h± and 
the component h z along the c axis are both nonzero. For 
a while we assume that the magnetic field h is given by 

h = h±e 2 + h z e 3 . (33) 

The explicit form of the potential of Eq. ((TJ) becomes 

2 , /a \2 



V 



1 r 
2 



+ 7r7 2 ™3 + o^i n 2 + h±h z n 2 n 3 + h±d z m , 



(dmY + (d 2 nY + lj - (34) 

[(dim - d 2 n 2 ) n 3 - (mdi - n 2 d 2 ) n 3 ] 

1 2 2 1 

— 7 ni H — 

2 ' 3 2 

where the parameter 7 2 depends upon /i z , 

7 2 = 1 + ^. (35) 

When h± 7^ 0, a brief inspection of the potential 
of Eq. |34|) reveals that the Zeeman energy | (n • h) 
now contains also the off-diagonal anisotropy h±h z n 2 ns, 
which was absent when either h± — or h z — 0. The 
presence of the latter anisotropy precludes analytical 
treatment. We therefore obtain the corresponding solu- 
tions by a direct minimization of the energy functional, 
in a manner analogous to the calculation presented in 
Sec. II. We state our results in the phase diagram in 

Fig.m 

We begin our discussion with the case where h z < 
h c i = 1.01 (or H z < H c i = 1.7 T). and consider the 
evolution of the system with increasing h± . Our results 
are displayed in Fig. HJa). In the limit hx = 0, shown 
in the top entry, the spin configuration that minimizes 
the energy is the flat spiral constructed in Sec. II. Recall 
that this solution is degenerate with respect to rotations 
around c, in agreement with the U(X) symmetry given by 
Eq. Q. When h± 7^ 0, our results make it clear that the 
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Figure 7. Theoretical prediction for the phase diagram. We adopted conventions used in publication of experiments 17 . The 
Antisymmetric phase is realized below the solid line. The Symmetric phase, denoted as S in the figure, exists in the area 
between the solid line and the dashed line. The dashed line depicts the limit of local stability of the Spin-flop phase. The 
Spin-flop phase is locally stable above the dashed line, but is actually realized only in the area above both the dashed line and 
the solid line. The inset emphasizes the region near the tricritical point V where the three phases (Symmetric, Antisymmetric, 
Spin-flop) merge. 



energy is minimized by a nonflat spin spiral propagating 
along the x axis. The component of the staggered mag- 
netization 72.2 is now different from zero and points in the 
direction of the transverse field h± , and its sign oscillates 
over the period L = L(h z ,h±). Because of this charac- 
teristic behavior, we call this state the Antisymmetric 
phase. The origin of the oscillating component ni can be 
understood by a direct inspection of the Zeeman energy 
oc (n-h) 2 . Its diagonal terms n\h\, n\h\ are always 
positive, but the off-diagonal contribution may become 
negative provided that 712 adjusts so that its sign is al- 
ways opposite to the sign of 773. Because the sign of 773 
oscillates (the projection of n onto the xz plane rotates 
during the period L), n-i also displays oscillatory behav- 
ior. To fully describe the spin structure, all terms in the 
potential of Eq. (fM|) must be considered, but the main 
conclusion persists - the spiral minimizes its energy by 
developing 772 ^ along the direction of the transverse 
field h±, and the sign of 772 oscillates over the period L. 
As a result, the expectation value (712/13) becomes nega- 
tive ( (712713) < 0), while (772) = (773) = 0, as verified by 
a direct calculation. 

We now briefly describe the role of the term h±d z n\ in 
Eq. (|34|) . The importance of the latter contribution has 
already been established in Sec. Ill during our analysis 
of the properties of the flat spiral (772 = 0) in the pres- 
ence of a field applied strictly in the xy plane (h± ^ 0, 
but h z = 0). The scenario discussed in Sec. Ill is here 
mildly modified by the presence of h z =/= but it main 
features remain the same, as confirmed by our numeri- 
cal studies. The weak-ferromagnetic anisotropy h±d z ni, 



generated by the transverse field hx applied along the y 
axis, makes the spin orientations along the ±x axis en- 
ergetically nonequivalent. In the Antisymmetric state, 
the component of the staggered magnetization that is 
perpendicular to the transverse field h± rotates in the 
xz plane, and is thus directly affected by the weak- 
ferromagnetic term h_\_d z n\. In turn, the profile of the 
Antisymmetric spiral is modified, and the expectation 
value of 77i over the period L becomes nonzero and neg- 
ative ((771) < 0) in order to minimize h±_d z n\. 

With increasing h±, the spiral becomes significantly 
distorted, and the 771 ~ — 1 orientation (domain) during 
the spin rotation is greatly enhanced. This is apparent 
from the bottom entry of Fig.|SJa). At the same time, the 
period L of the spiral increases, and the energy density of 
the Antisymmetric state begins to approach the energy 
density of the uniform Spin-flop state n = (—1, 0, 0) from 
below. At the critical value of the transverse field /i]_(/i z ), 
the period of the spiral grows to infinity [L — > 00), and 
its energy density becomes equal to the energy density 
of the Spin-flop state w — h (1 — 2h±d z ). This numeri- 
cally verified scenario is consistent with the experiment^, 
and is somewhat similar to that discussed in Sec. Ill for 
strictly transverse fields. Above the critical line, only the 
uniform Spin-flop states emerges from our numerical cal- 
culations, and the incommensurate Antisymmetric spiral 
no longer exists. The boundary between the Antisym- 
metric and the Spin-flop state is indicated by the solid 
line in Fig. [7] We have verified that the Antisymmetric 
state displayed in the phase diagram always carries lower 
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(a) h z =0.59 (b) h ± = 2.68 




Figure 8. Calculated evolution of spin configuration in the 
Antisymmetric phase with the applied magnetic field, (a) Left 
panel: h z = 0.59 (H z = 1 T), h± increases, (b) Right panel: 
h± = 2.68 (H± = 4.5 T), h z increases. The bottom entries 
in both panels are applied for points near the critical line 
between the Antisymmetric and the Spin-flop phase. Notice 
an enhanced m = — 1 domain in these entries. 



energy density than the uniform Spin-flop state. 

Evolution of the spin structure with increasing h z , but 
fixed strength of the transverse field h±, is shown in 
Fig.EJb). Our results, applied here for h± = 2.68 (H± = 
4.5 T), are qualitatively similar to those in Fig. E[a). In 
particular, the spiral again develops a nonzero oscillat- 
ing component ri2 7^ along h±, with zero expectation 
value (712) = over the period L. Expectation value 
(^2^3) < due to the off-diagonal anisotropy n2n^h±h z , 
whereas (n\) < thanks to the weak-ferromagnetic term 
h±d z ni. With increasing h Zl the period of the spiral L 
increases and presumably again diverges (L — » 00) at 
the critical line. Above the critical line, only the uniform 
Spin-flop states emerges from our numerical calculations, 
and the incommensurate Antisymmetric spiral no longer 
exists. 

We emphasize, that the characteristic properties of 
the Antisymmetric state discussed in the preceding para- 
graphs remain the same for any point h z ^ 0, hj_ ^ 
below the solid line in Fig. [7] However, the scenario 
of the phase transition between the Antisymmetric and 
the Spin-flop phase, discussed in connection with Fig. HI 
is slightly modified for sufficiently weak h±, near the 
point r. Specifically, for H± below ~ 1 T, the energies 
of both states again become equal at the critical line, 



but the period L of the Antisymmetric spiral remains 
finite (albeit large). Above the critical line, our numer- 
ical minimization still yields the solution in the form of 
the Antisymmetric spiral, with however, the energy den- 
sity higher than the energy density of the Spin-flop state 
n = (—1,0,0). This should be contrasted with behavior 
for large h±, where the period of the Antisymmetric spi- 
ral diverges (L — > 00) at the critical line; and above the 
critical line only the uniform Spin-flop state exists. In- 
terestingly, our results seem to be again consistent with 
the experiment^. 

To complete our description of the Antisymmetric 
state, we display the typical 3D path traced out by the 
endpoint of spin during one period L in Fig.lTOfc), calcu- 
lated for h z = 1.21 and h± = 0.28/1* = 0.34 (H z = 2.04 
T and H± = 0.57 T). The spin approximately rotates in 
a plane whose normal is tilted away from the y axis in 
the yz plane. 

To proceed further, we now explain the meaning of 
the dashed line in Fig. [71 and examine an important 
issue concerning the existence of the Spin-flop state in 
the presence of arbitrary canted fields. We note that 
the uniform Spin-flop state n = (—1,0,0), or $ = —5, 
Q = ^ using the spherical parametrization ([S]), is an 
essentially obvious stationary point that minimizes the 
energy functional W = J Vdx dy, where V is the poten- 
tial given in Eq. (f3~4"|) . Actually, there exist two differ- 
ent spin-flop configurations n = (=f1,0,0) and both of 
them are the stationary points of the corresponding en- 
ergy functional. However, their energy densities given by 
W = I (1 T 2h±d z ) are different. Therefore, we will only 
consider the Spin- flop state n = (—1, 0, 0) with lower en- 
ergy in our analysis. To examine the stability, we first 
introduce new fields 

y,t) = -- + f{x, y, t) , 6(2;, y, t) = - + g(x, y, t) , 

(36) 

where f(x,y,t), g(x,y,t) account for small fluctuations 
around the Spin-flop state. Now the actual parametriza- 
tion of the staggered magnetization n given by Eq. (|3"6"]l 
is inserted in the complete Lagrangian of Eq. ((T|), which 
is applied for a magnetic field h given by Eq. (|3"3")l and ex- 
panded to quadratic order in /, g. If we further perform 
the usual Fourier transformation with frequency u and 
wave vector q = (91,(72), the corresponding linearized 
equations of motion can be solved analytically to yield 
the (squared) eigenfrequencies 

u4(q) = $ + q$ + h±d z + (37) 

+ \ + h 3 z + h\ ± ^{l + hl-h\f +4/^2+16^ . 

The above calculated magnon spectrum is strongly 
anisotropic. To examine the local stability of the Spin- 
flop state, we note that the stability condition requires 
that uPl > and uPl > for each q. It is also clear that 
uj\ > and uPl is minimum for qi — 0. Therefore, we 
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0) 



h x /h z 



0.08 



(b) 



h ± /h z 



0.28 



3 -hi -2\h 



4h±d z . 



(38) 



The above obtained line of local stability of the Spin-flop 
state is displayed by the dashed line in the phase 
diagram of Fig. [7J Below the dashed line, the Spin-flop 
state is locally unstable. Thus the Spin-flop state cannot 
exist beyond the point T, where the energy density of 
the Antisymmetric state becomes equal to the energy 
density of the Spin-flop state. 

It is more or less clear, there is a new phase realized 
in some area just below the dashed line, near the axis 
H z (near H± =0). Note that Eq. ([38]) applied for the 
special case h± = yields h z = -J3 (H z = 2.9 T), which 
is just the upper critical field h C 2 obtained in Ref. 12. For 
h z < h C 2, the Spin-flop phase is locally unstable, and the 
intermediate phase is realized in the region h c \ < h z < 
h C 2- It is then natural to expect that the intermediate 
phase exist in some form also in the presence of a weak 
transverse field h± 7^ 0. 

Our calculations confirms this expectation. Specif- 
ically, when h± 7^ 0, the so-called Symmetric phase 
emerges as the ground state in the region between the 
dashed and the solid line in Fig. [7J Examples of the 
Symmetric spiral, calculated for two values of the ap- 
plied magnetic field, are illustrated in the middle and 
the bottom entry of Fig. [9ja). The typical path, traced 
out by the endpoint of the spin is shown in Fig. [TWb'l . 
The Symmetric phase can be described as an antifcrro- 
magnetic conical spiral that propagates strictly along y, 
but nutates around the —x axis. Importantly, the com- 
ponent ni, perpendicular to the transverse field hj_, is 
nonzero (nj 7^ 0) and always negative ((rii) < 0). All 
these features agree with the experiment^. 

The Symmetric phase is a direct descendant of the in- 
termediate phase - the conical antifcrromagnetic spiral- 
discussed in Sec. II for fields strictly parallel to c (h± = 
0). Recall that the intermediate phase obeys the U(l) 
symmetry described by Eq. ©. In practice, the U(l) 
symmetry is broken by an additional tetragonal sym- 
metry induced by discreteness effect o 6 ' 1. Thus, in the 
absence of transverse fields, there exist four degenerate 
states, illustrated in Fig. [TOTa): the intermediate spiral 
propagates along x and nutates around the ±y axis; or 
it propagates along y but nutates around the =px axis. 

The above described degeneracy is broken when h± 7^ 
0. To see that, we consider for a moment the interme- 
diate spiral n calculated for h± = 0, and insert this so- 
lution in the potential (IM)) applied for h± 7^ 0. The 
additional corrections to the energy are then given by 
g^x^i) + h±d z (ni). The first correction, quadratic in 
hj_, originates in the Zeeman energy oc (n-h) . Note 
that the off-diagonal term h±h z (n2n^) in the Zeeman en- 
ergy does not contribute, because the expectation value 
(712/13) vanishes in the intermediate phase for any degen- 
erate state. The second correction, linear in h±, is due 




~\ h z =0.70 


/ 3 A 








Figure 9. Examples of spin configurations in various 
phases from the phase diagram in Fig. [7J calculated for (a) 
h±/h z =0.08, and (b) h±/h z = 0.28. The two values h±/h s 
correspond to the magnetic field applied at an angle 4.57° and 
15.64° with respect to the c axis. The middle and the bottom 
entry in (a) corresponds to the Symmetric phase, where the 
spiral propagates along y. All other entries show the Antisym- 
metric phase, with propagation direction along x. Phase tran- 
sition between the Antisymmetric and the Symmetric phase 
is accompanied by sudden 7r/2 rotation of propagation direc- 
tion. 



to the weak-ferromagnetic term d z (h x e^) ■ n. The lin- 
ear weak-ferromagnetic contribution dominates for small 
transverse field, and favors the conical spiral propagating 
along y and nutating around the —x axis, with (ni) < 1. 

Numerical work confirms that the above argument is 
correct despite the simplifying assumption that neglects 
the changes in the staggered magnetization induced by 
h± =/= 0. Thus, the Symmetric phase illustrated in 
Fig. [91(a) and Fig. [TOlbb is essentially the intermediate 
phase of Sec. II, propagating along y but nutating around 
the —x axis, whose profile n is modified by h± 7^ 0. 



Sudden rotations 

The Symmetric phase emerges in canted magnetic 
fields applied nearly parallel to the c axis, when h z > h c \. 
It is the stationary point of the energy functional with 
the lowest energy density in the area between the dashed 
line and the solid line of Fig. [7J With increasing h z the 
magnitude of n\ becomes larger and larger, until at the 
dashed line n± — ¥ — 1 and the solution becomes the Spin- 



h z =1.21 

(a) h ± =0, Intermediate phase 

propagates along x, nutates around +y or -y 
z z 





propagates along y, nutates around -x or +x 
z z 





(b) h ± =0.1 , Symmetric phase (c) h ± =0. 34, Antisymmetric phase 



propagates along y 



propagates along x 





Figure 10. (Color online) Typical spin configurations in 
canted magnetic fields, calculated for h z — 1.21 and the three 
values of the transverse field h±. Thick lines on the sphere 
indicate paths traced out by the endpoints of the staggered 
magnetization during one period L. The base of the stag- 
gered magnetization is placed at the center of the sphere, (a) 
Intermediate phase, h± = 0. The conical spiral may propa- 
gate along x and nutate around the ±y axis; or may propa- 
gate along y and nutate around the =px axis, (b) Symmetric 
phase, h± =0.1. The conical spiral propagates strictly along 
y but nutates around the —x axis, (c) Antisymmetric phase, 
h± — 0.34. Spin approximately rotates in a plane, whose nor- 
mal is tilted away from the axis y in the 1/2-plane. The spiral 
propagates strictly along the x axis. 



flop state n = (—1,0,0). These results generally agree 
with experimental finding o 16 ! 17 . Evolution of the spin 
structure with increasing hj_, but fixed strength of the 
longitudinal component h z is rather mild. The period 
L slightly decreases with h±, whereas the magnitude of 
ni moderately increases due to the weak-ferromagnetic 
energy d z (h x e 3 ) ■ n. 

Importantly, at the critical solid line the energy den- 



B 
"S 

E 

CO 



1.0 



0.8 



CO 
CL 

If 0.6 



3 4 

CD 

1 0.2 
o 
o 
c 

0.0, 



T TT 77 1 



~i — 1 — i — 1 — i — 1 — r 



• o 



o H ± /H z =0.08 
. H ± /h)=0.28 



j i i i i i i_ 



0.5 1 1.5 2 2.5 3 3.5 4 
H Z (T) 

Figure 11. Theoretical field dependence of the incommensu- 
rability parameter calculated for H±/H z =0.08 and 0.28. The 
two values H±/H z correspond to the magnetic field applied 
at an angle 4.57° and 15.64° with respect to the z-axis. Our 
choice closely follows values 5° and 15° used in the experi- 
mental Fig. 10, Ref. 17. The SI units are used to facilitate 
comparison with the experiment. For H±/ H z =0.08, a pre- 
dicted discontinuous jump of the incommensurability param- 
eter by ~ 5% at H z w 1.89 T (dashed line) corresponds to 
the phase transition between the Antisymmetric and the Sym- 
metric phase. This is consistent with the experimental data. 
The dotted line corresponds to the critical field H c i. 



sity of the Symmetric phase becomes equal to the en- 
ergy density of the Antisymmetric phase. For the spe- 
cial choice h z — 1.21 shown in Fig. [TUl this happens at 
h± m 0.13 or 0.22 T. For stronger h±, the Antisymmetric 
state (Fig. fTPT c)) emerges as the true ground state, with 
the energy density lower than the Symmetric spiral . 

Because the Symmetric spiral propagates strictly along 
y, whereas the Antisymmetric phase propagates along 
the x axis, the above described first order phase tran- 
sition is accompanied by a sudden rotation of the spiral 
propagation direction exactly by 5. This theoretical re- 
sult agrees with experiment o 16 i 17 . 

In order to compare with the experiment directly, we 
present our results calculated for two values of the pa- 
rameter h ± /h z = 0.08 and/or 0.28, in Figs. H [Til and 
[T2"l Our choice of h±/h z corresponds to a magnetic field 
applied at the angle a = 4.57° and 15.64° with respect 
to the c axis, which is similar to the angles 5° and 15° 
used in neutron scattering measurements 17 . 

Fig. HJa) applied for fields nearly parallel to the c axis 
(h±/h z — 0.08, or a = 4.57°), summarizes the evolution 
of the spin structure with increasing strength of the field 
across the phase transition. For h z = 0.70 < h c \ (the top 
entry), the system is in the Antisymmetric phase. The 
Antisymmetric spiral propagates along x, perpendicular 
to the transverse field h± . Its oscillating component n2 is 
for given value of the field rather small, and the structure 
resembles the flat spiral (712 = 0). It may be identified 
with the soliton lattice of Ref. 17. For h z > h c \, the 
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Figure 12. Theoretical field dependence for the intensities of 
the 1st, 2nd and the 3rd Fourier harmonics calculated from 
the ni and 713 components of the staggered magnetization. 
We used H±/H z =0.28. The latter parameter corresponds to 
the magnetic field applied at an angle 15.64° with respect to 
the z-axis and is roughly equal to 15° used in the experimental 
Fig. 8, Ref. 17. We adopt SI units to facilitate comparison 
with the experiment. 



Symmetric state - a conical structure propagating along 
y, but nutating around the —x axis emerges (the middle 
and the bottom entry). Note that h± is applied along y, 
and our result is thus in agreement with Fig. 13 of Ref. 17. 
The Symmetric spiral should be identified with the AF 
cone of Ref. 17. The phase transition at the critical field 
h z w h c i is of first order, and is accompanied by the 



sudden J rotation of the spiral propagation direction. 
This experimental observatio n 16 ' 17 is explained by our 
theoretical analysis. The magnitude of ni increases with 
the strength of the field, and m — > — 1 at the critical field 
w 2.7 T. For stronger fields, the system is in the Spin-flop 
state. The calculated value of the IC phase transition is 
slightly larger than the experimental value 2.4 T. 

The calculated field dependence of the incommensura- 
bility parameter L(0)/L(h) is presented in Fig. [TT] (empty 
diamonds). Note a discontinuous jump at ~ 1.89 T where 
the incommensurability parameter increases by rs 5%. 
This corresponds to the phase transition from the An- 
tisymmetric to the Symmetric phase and highlights its 
first order nature. We emphasize a remarkable agree- 
ment with experimental data in Fig. 10 of Ref. 17 and 
quoted values 1.95 T and ~ 6% for the critical field and 
a discontinuous jump of the incommensurability param- 
eter. Actually our data are shown only for field values 
up to 2.6 T < H C 2 ~ 2.7 T. However, our calculations 
indicate a continuous IC transition with n\ — > — 1, but 
finite L in the limit H —¥ H c i- 

We now discuss the case where the field is applied at 
"large" angle a s» 15° (h±/h z = 0.28) with respect to 
the c axis. Evolution of the spin structure with increas- 
ing strength of the field is summarized in Fig.UJJb), while 
the corresponding field dependence of the incommensu- 
rability parameter L(Q)/L(h) is illustrated in Fig. HT1 (full 
circles). In this case, no reorientation of the spiral propa- 
gation direction was observed in the experiment^ 7 -. This 
is consistent with our T = calculations, which predict 
the Antisymmetric phase with oscillating 712 (paral- 
lel to h±) propagating along x for all field strengths until 
the IC transition at ps 3.33 T. Predicted critical field is 
somewhat larger than observed in the experiment, but is 
not terribly inconsistent with the value 2.6 T quoted in 
Ref. 17. The 712 component is small for weak fields, but 
its magnitude quickly increases with the strength of the 
field. For strong fields, the structure becomes clearly non- 
sinusoidal, and non-planar and can be identified with the 
"distorted incommensurate structure" of Ref. 17. This is 
seen in the middle and the bottom entry in Fig. [9jb). 
The Fourier transform of the staggered magnetization 
provides evidence for higher harmonics, both odd and 
even. 

Our T = theoretical results for the field dependence 
of the intensities of the 1st, 2nd and 3rd Fourier com- 
ponents of the staggered magnetization n(x) are pre- 
sented in Fig. [T3J The intensities are calculated from 
the Tii, 773 components. This is because neutron scatter- 
ing sees only the components perpendicular to momen- 
tum transfer (which is parallel to 712). Our results are 
to be compared with experimental data in Fig. 8(b) of 
Ref. 17. For the 1st harmonic the agreement with ex- 
periment is apparent, except that we predict somewhat 
higher value for the of the IC critical field. Higher har- 
monics smoothly appear above 1.7 T - 2 T, and the in- 
tensity of 2nd harmonics linearly increases with the field, 
as in the experiment. Similarly, the 3rd harmonics first 
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increases, then shows a shallow dip and increases again 
close to the IC critical field. This is again fully consistent 
with experimental observation. On the other hand, our 
results show rapid increase of both higher harmonics as 
the field approaches the critical value, whereas the exper- 
imental data show smoothing at the IC transition. This 
can be perhaps due to finite temperature. Overall agree- 
ment with experiment is however fairly good. Finally we 
discuss the field dependence of the incommensurability 
parameter L(0)/L(h) shown in Fig. [TT] (full circles). In 
agreement with experiment, no discontinuous jump in the 
theoretical prediction is observed. The curve is concave 
for weak field strengths, but becomes convex for if > 1.7 
T. This corresponds to the emergence of higher harmon- 
ics and is again in agreement with the experimental data 
in Fig. 10 of Ref. 17, with minor discrepancy in the value 
of the IC critical field. 

We end this section with two comments: 

• Numerical work confirms that the existence of d z ^ 
is not crucial for the appearance of the double- 
k structure and/or sudden 7r/2 rotations observed 
in experiment. It is, however, important to restore 
quantitative agreement with experiment. In partic- 
ular, for d z — and fields nearly parallel to c, the 
intermediate phase would first appear at h c ± in the 
form of a nonflat spiral propagating strictly along 
x but nutating around the y axis, without reorien- 
tation of the spin propagation direction. A sudden 
7r/2 rotation of the propagation direction occurs 
later, at yet another critical field h ro t a tion ~ 1-20, 
above which the minimum energy state becomes 
the Symmetric spiral propagating along y but nu- 
tating around the x axis. In the absence of the 
weak-ferromagnetic energy d z (h x e^) ■ n, an expla- 
nation of sudden reorientation requires a detailed 
analysis of the of the energy term oc (n • h) 2 . 

• We assumed that the transverse component of the 
field h±_ points strictly along the y axis. Our re- 
sults, however, are not restricted to this special 
case. For example, assume that hj_ points in an 
arbitrary direction in the xy plane, which is ob- 
tained by a clockwise rotation of the y axis with 
angle Then the staggered magnetization n for 
any state calculated earlier in this section must be 
also rotated clockwise with the angle ip around the 
c axis, while the original direction of spin propa- 
gation must be rotated counter-clockwise, with the 
angle —ip. All other results remain unchanged. 



V. CONCLUSION 

We have presented a rather complete theoretical study 
of T = phase transitions in canted fields of arbitrary 
strength and direction. We calculated the complete phase 
diagram and identified the symmetries of states in a num- 
ber of different regions. For the fields applied nearly 
parallel to the c axis, we confirmed the existence and 
stability of the "intermediate" phase that mediates the 
incommensurate-commensurc transition and analyzed its 
properties. We identify this phase with an experimen- 
tally observed double-k structure. By analyzing data on 
fields applied perpendicular to the c axis, we determine 
an out-of-plane anisotropy parameter d z needed to com- 
plete quantitative comparison with experiment. Finally, 
our model accounts for sudden 7r/2 rotations that have 
been highlighted as a noteworthy feature of recent exper- 
iments. 

The work reported in this paper results from a 
long-standing theoretical investigation of spiral mag- 
netic structures in Dzyaloshinskii-Moryia antiferromag- 
nets. The theoretical framework involves a number of 
approximations: the replacement of quantum-mechanical 
by classical variables, ignoring inter-layer couplings and 
the replacement of discrete spins by continuous fields in 
a model Lagrangean. 

Nevertheless, detailed agreement with experimen t 16 ! 17 
is now so extensive that the applicability of this model to 
systems such as Ba2CuGe20y may now be established. 
The only remaining discrepancies lie in the particular 
magnetic field values at which transitions between mag- 
netic states take place, and these discrepancies are on the 
order of 10-20%, which is not much beyond experimental 
uncertainty. 
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